Inconsistency in the application of the adiabatic theorem 
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The adiabatic theorem states that an initial eigenstate of a slowly varying Hamiltonian remains 
close to an instantaneous eigenstate of the Hamiltonian at a later time. We show that a perfunctory 
application of this statement is problematic if the change in eigenstate is significant, regardless of 
how closely the evolution satisfies the requirements of the adiabatic theorem. We also introduce an 
example of a two-level system with an exactly solvable evolution to demonstrate the inapplicability 
of the adiabatic approximation for a particular slowly varying Hamiltonian. 

PACS numbers: 03.65.-w, 03.65.Ca, 03.65.Ta 



Introduction. — Since the dawn of quantum mechanics 
miUlii' t ne venerable adiabatic theorem (AT) has under- 
pinned research into quantum systems with adiabatically 
(i.e. slowly) evolving parameters, and has applications 
beyond quantum physics, for example to electromagnetic 
fields. The AT lays the foundation for the Landau-Zener 
transition (LZT) (including the theory of energy level 
crossings in molecules) |4| , for the Gell-Mann-Low theo- 
rem in quantum field theory || on which perturbative 
field theory is constructed, and for Berry's phase Q. 
More recently the AT has renewed its importance in the 
context of quantum control, for example concerning adia- 
batic passage between atomic energy levels , as well as 
for adiabatic quantum computation ||. Essentially the 
standard statement of the AT implies that a system pre- 
pared in an instantanous eigenstate of a time-dependent 
Hamiltonian will remain close to an instantaneous eigen- 
state of the Hamiltonian provided that the Hamiltonian 
changes sufficiently slowly. Here we demonstrate that the 
application of the standard statement of the AT leads to 
an inconsistency, regardless of how slowly the Hamilto- 
nian changes; although the AT itself is sound provided 
that deviations from adiabaticity are properly accounted 
for, the standard statement alone does not ensure that 
a formal application of it results in correct results. In 
addition, we will present a simple two-level example for 
the failure of the AT even in a case when all the criteria 
for the AT seem to be met. The results of this paper are 
intended to serve as a warning that the incautious use 
of the AT may produce seemingly unproblematic results 
which nevertheless may be grossly wrong. 

Before demonstrating the inconsistency arising from 
the standard statement of the AT, it is useful to give 
a simple exposition of the theorem's proof as given in 
Ref. [2j . The evolution of a quantum state \4>(t)) un- 
der a unitary evolution U(t,to) is described by \t/)(t)) — 
U(t,to)\Eo(to)). We assume that the initial state 
\Eo(to)) is an eigenstate of the initial Hamiltonian H(to). 
The time-dependent Hamiltonian is related to U(t) by 
H(t) = WW 18], such that |V>(*)> fulfills the usual 
Schrodinger equation. In the instantaneous eigenbasis 



{\E n (t))} of H(t), the state can be expressed as \i/j(t)) = 
Yl n ' t l ; n{t)e~ 1 I E "\E n (t)), where we have introduced the 
short-hand notation J E n = J t E n (t')dt' . Inserting this 
expansion into the Schrodinger equation leads to the fol- 
lowing differential equation for the coefficients, 



iip n = -i ^2 e ' 



f(E»-E m ) 



■tp m {E n \E m ) 



(1) 



The AT relies on the requirement that H(t) is slowly 
varying according to 



\{E n \E m )\ « \E n -E„ 



n ^= m 



(2) 



Transitions to other levels are then supposed to be neg- 
ligible due to the rapid oscillation arising from the phase 
factor exp(i J (E n — E m )), yielding 



W)) 



-IS Bo. 



[0 °\E o (t)) 



(3) 



with p n = i J (E n \E n ) the geometric phase (GP) jy]. 
Condition J5J and approximation © summarize the 
standard statements of the AT. 

Proof of inconsistency. — The inconsistency implied 
by Eq. (0) is evident by considering the state := 
UHt,t )\E (t )). Using d t (UW) = it is easy to see 
that this state fulfills an exact Schrodinger equation with 
Hamiltonian H(t) = -W (t,t )H(t)U(t,t Q ). To demon- 
strate the inconsistency, we commence with a claim that 
is shown to yield a contradiction. 
Claim: The AT © implies 



$} = e i I E °\E (t )) 



(4) 



Proof of inconsistency: Because U(to,to) = 1, re- 
sult 0} fulfills the correct initial condition so it remains 
to show that also fulfills the Schrodinger equation: 

idS) = -E (tM 

= -E {t)U^Ue^ Eo \E {t )) 

w -E Q {t)U^e^\E {t)) 

= -U^H(t)e if3 «\E Q (t)) 

» -tfH{t)Uj$ E °\E {t )) 

= mm ■ (5) 
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The AT is explicitly used in the lines with 
Eq. J3J implies 



However, 



(E Q (t )\Utf\E (t )) = (Eo(to)\U$) 

« e tA ><JS^(* )|-Eb(*)> ^ 1 , (6) 

which is false □. 

Clearly the inconsistency is a consequence of neglecting 
the deviations of Eq. © from the exact time evolution 
which is free of inconsistencies. Stated another way, ap- 
proximation Pjl. without correction terms, could only be 
exact in the limit of infinitesimally slow evolution, for 
which the system is constant over finite time and the 
evolution is indeed given by a multiplicative phase cofac- 
tor. However, evolution is not infinitesimally slow, and 
neglect of the correction terms leads to the inconsistency 
demonstrated above. To elucidate this point we define 
the following unitary transformation, 



U A T(t,t ) 



-'tto^e^lE^tyiEMl . (7) 



The (exact) time evolution generated by Uat is equiva- 
lent to the standard statement © of the AT for adiabatic 
motion in a finite-dimensional Hilbert space with non- 
degenerate energy levels. It is straightforward to write 
-Hat(£) = — iU AT U at in the form 

#at(*) = -^2E n (t)\E n (t ))(E n (to)\ 



,iS{E n -E m )-i{0 n -l3 m ) 



x{E n {t)\E m {t)) \E n (t ))(E m (t )\. (8) 

The second sum in this expression has the same structure 
as those terms that are omitted in the adiabatic approxi- 
mation, and by omitting these terms one again arrives at 
the inconsistent result Q. However, evaluating Hamil- 
tonian -Hat hi the interaction picture with respect to the 
first line of Hat, one finds 

H AT (t) = -i J2 e-^-M 

x(E n (t)\E m (t)) \E n (t ))(E m (t )\ . (9) 

In this Hamiltonian the transition matrix elements be- 
tween the different initial eigenstates are not rapidly 
oscillating anymore and therefore cannot be neglected. 
However, perfunctory use of the standard statement of 
the AT © implicitly neglects such terms. 

Thus we have shown that the standard statement of 
the AT may lead to an inconsistency no matter how 
slowly the Hamiltonian is varied, but so much science 
rests on the AT that the implications of this inconsis- 
tency are important and require exploration. Perhaps 
the most important application of the AT is the slow evo- 
lution of an initial instantaneous eigenstate \Eo(t)) into 



a later instantaneous eigenstate \Eo(t)) that is meant to 
be quite different; i.e. = \(E {t)\E (t ))\ < 1. For 
example the famous LZT evolves a two-level molecule 
or atom with orthogonal basis states |0) and |1) from 
\Eo(to)) = 10} to \E {t)) — |1) with near-unit probability 
so that !Fq « |(0|1)| = 0. On the other hand, the quan- 
tity Tx = \(E a (0)\UW\E o {0))\ should always be unity, 
but © implies T\ ~ « 0. Thus, the deviation of the 
overlap function T§ from unity is an alarm indicator for 
when the AT is vulnerable to the inconsistency: whenever 
\E (t)) deviates strongly from the initial state |i?o(0)}, 
the inconsistency is a potential problem, regardless of 
how slowly H(t) changes. 

Counterexample of a two-level system. — The incon- 
sistency introduced above is due to a particular inverse 
time evolution which causes rapidly oscillating terms to 
become slowly varying. One may see this as a resonance 
problem which can also appear for U itself. As a specific 
example, consider a two-level system with exact time evo- 
lution defined by 

U(t) = exp (-i0(t)n(t) • er) = cos 91 -in -a sin 9 (10) 

with 6(t) = ujot, n(t) = (cos(27rf/r), sin(27rf/r), 0) and 
er = (cr x ,&y,&z) denoting the Pauli spin vector oper- 
ator. The associated Hamiltonian can be calculated 
using H(t) — iUU^ and can be written in the form 
H(t) = R(t) • er, with 



R 



9n + cos#sin(9ri + sin 2 6(n x h) 



LL> Q n(t) 



27rsin(ti;o^) 



_ sin (2^)cos(w i) \ 
C08(^*) COS((J t) 

sm(ujQt) J 



, , sin(cjni) -a / x 
uj n{t) + — K -^K(t) 



(11) 



This Hamiltonian is similar to that of a spin-i system in a 
magnetic field of strength proportional to loq that rotates 
with period r in the x — y-plane. The exact Hamiltonian 
for the latter case would correspond to R(i) = 0; we 
will discuss the importance of this difference below. The 
eigenvalues of H (t) are given by 



E±(t) = ±|R(f)| = ±\J9 2 + sin 2 On 



(12) 



It is easy to show that the evolution operator (|1U|) fulfills 
requirement (0) for adiabatic evolution as long as the 
vector n changes slowly compared to ujq, i.e. for u t ^> 1. 
The time scale r corresponds to the large time scale which 
appears in the mathematically more elaborated forms of 
the AT 0, 0, 0| ■ These correction terms are resonant 
[l^| so that a large deviation from the AT predictions can 
accumulate over time. 

To evaluate the predictions of the AT it is convenient 
to consider projection operators instead of the state itself. 
Projectors onto eigenstates of H(t), which fulfill H(t) \ ± 
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(t)) = ±|R(t)| | ± (t)), can generally be written as 



P 



l±(t)> 



R(t)| 



(13) 



If we consider the evolution at time T — t/2 and as- 
sume for simplicity that ojqT is a multiple of 2tt, we have 
R(T) = -R(0) and C/(T) = 1. We thus find P|+ ( t)> = 



|-(0)> 



but iV|+(o)) - t/(T)P| +(0)) C/(T)t = Pj. 



(o)>- 



In 



other words, the perfunctory prediction 

U(T)P l+m U(T)i « P, +(r)) 



(14) 



of the AT is invalid. Thus, whereas a resonant but weak 
time-dependent oscillatory term in the evolution repre- 
sents an unusual application of the AT, this system meets 
the criteria of the AT and therefore casts doubt on the 
general applicability of criterion (0). 

For two- level systems, it is possible to derive a general 
criterion on when the AT is bound to fail, i.e., when the 
quantity Q := \(+(t)\U(t)\ + (0)}| 2 - Tr P c/ | +(0)) P| +(i)) 
strongly deviates from one. It is evident that this crite- 
rion depends on U (t) at time t only, as well as on the 
Hamiltonians H(t) and H(Q). There is no direct refer- 
ence to the slow evolution of the Hamiltonian because 
the criterion does not depend on H. For a unitary trans- 
formation of the form I jlOjl with general 9(t) and n(i) it 
is straightforward to derive 



l + n(0) 



9n + cos 9 sin 9h — sin 2 On x h 

iRi 



(15) 

We have assumed that U(0) is given by the identity ma- 
trix so that 9(0) = and R(0) = #(0)n(0). To examine 
when Q can become small we focus on a special case 
of adiabatic evolutions, characterized by 9 3> |n|. In 
this case we can neglect all terms containing ri such that 
R| pa 8. We then arrive at the conclusion that the AT 
is maximally violated if n(t) w — n(0), as in the case for 
the example given above. We remark that many other 
adiabatic evolutions do not fulfill 9 » |h|, since it implies 
that 9 > so that U (t) has to become equal to the iden- 
tity again within the fast time scale 1/9. For instance, a 
LZT, for which R(t) = ile x — A t/2e z for constant real 
f2 and A , asymptotically fulfills R(t) ~ — R(— t), but 
not 0> |h| so that Q « 1 is still valid. 

Eq. (|15f) is a universal criterion for the failure of the 
AT for two-level systems. Although it is likely that the 
small but resonant terms in our counterexample (CE) 
are the cause for this failure, a non-resonant CE to the 
consistency of the standard statement of the AT is not 
necessarily excluded. This is because Eq. (|15|l depends 
only on the initial and final Hamiltonian and the final 
unitary matrix U, and therefore makes no reference to 
the behaviour during the evolution. 

It is worthwhile to examine if a resonant behaviour as 
in our CE is excluded by the conditions imposed on the 



Hamiltonian in more rigorous forms of the AT. The two 
cases we do consider both require the usual gap condition 
for the energy levels, which is fulfilled in the CE. In addi- 
tion, Kato [3j demands dH(s)/ds to be finite for r — » oo, 
where s = t/r is a scaled time variable. This is the case 
for the CE |2fJ. In another proof of the AT, Avron et 
al. |l(j require the Hamiltonian to be at least twice con- 
tinuously differentiable, which is also fulfilled by the CE 
[1H . In this case the AT (Theorem 2.8 of Reference H3) 
is slightly different and states that Pu\+(o)) stays close to 
P;7 A |+(o)), where the unitary operator UA(t) is generated 
by the modified Hamiltonian 



H A (t) = H(t) + i[P ] 



+(*)>> -H+(*)>J 



(16) 



(c.f. Eq. (1.0) and Lemma 2.2 of Ref. 10]). For the 
CE presented above, we have numerically solved the 
Schrodinger equation (in the scaled time s = t/r) for the 
propagator Ua and calculated the fidelity (or overlap) 



T = TrJP, 



,1/2 



u\+(o)) Pu A \+(o)) ^u\+(o)) 



p, 



1/2 



(17) 



between the exact time evolution and the eigenvector 
subspace propagated with Ha- The result is shown in 
Fig. As in our analytical results the overlap becomes 
zero for t/r = 1/2 where the maximal violation occurs. 

Thus it seems that the conditions on the AT are not 
strict enough to exclude the CE. A way to exclude reso- 
nant but small behaviour may be to demand continuous 
differentiability of H (s) even in the limit r — > oo. How- 
ever, while this would be a sufficient criterion to exclude 
resonances it may not be a necessary criterion and thus 
could exclude other cases in which the AT works well. 
Also, since it is not proven that resonances are the cause 
of problems, this criterion might not exclude other cases 
where the AT may become problematic. 

Remarks on the validity of the AT. — Although the 
standard statement of the AT may be problematic in cer- 
tain applications, previous results based on the AT are 
generally not necessarily affected. The reason is that the 
inconsistency is not related to the validity of the AT as 
an approximation but to its application in formal deriva- 
tions. In addition, most applications of the AT as an 
approximation do not include resonant perturbations, so 
that the AT should provide an excellent approximation 
to the exact time evolution. This is the case, for in- 
stance, for a real spin-^ system in a slowly rotating mag- 
netic field (R = in the Hamiltonian above) and for 
LZTs. The correctness of the LZT may also guarantee 
that the results of adiabatic quantum computation || re- 
main valid because, for a two-level system, the latter can 
be mapped to the first. However, if the reversed time 
evolution U^(t 7 t ) were to be computed using Eq. Q), 
the inconsistency could yield an incorrect state. 

An example where the inconsistency associated with 
the AT poses a significant problem is a perturbative 
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treatment of the GP. For brevity we refer to Refs. |g, 
EL 111) for explanations of the technical terms and the GP 
used in this paragraph. Under the condition of parallel 
transport the GP of an evolving state is given by the 
phase of (ip(0)\ip(t)) = (ip(0)\U(t)\ip(0)). If we consider 
the case that the unitary operator is slightly perturbed 
by an operator P, one can show that for an open quan- 
tum system the associated corrections include terms of 
the form (ip(0)\U(t)P\i()(0)) and <^(0)|Pl7(t)|^(0)> [H- 
In order to calculate these corrections one needs in par- 
ticular to find an expression for the state {tp(0)\U(t) = 
(W (t)\ijj(0))y . It is obvious that the inconsistency would 
then lead to a wrong result for the GP. 

In general, a potential problem in the application of 
the AT could be the presence of small fluctuations in an 
experiment, even if the ideal case would not be affected 
by the inconsistency. The reason is that example 1101) in- 
dicates that small changes can invalidate the predictions 
of the AT, even if they respect the adiabadicity criterion 
(J2J. In the two- level CE, the omission of the small terms 
proportional to R in the Hamiltonian changes a system 
where the AT is valid to one where it is maximally vi- 
olated. Likewise, the Hamiltonian JHJ shows that it is 
exactly the omission of the small terms which leads to 
the inconsistency. Thus whenever adiabatic fluctuations 
are present in an experiment, it seems to be necessary 
to check the predictions of the AT. This could be done 
by checking the quantities !Fq and Q for mixed states. 
To be more specific, we consider a system with fluc- 
tuations in the classical parameters that determine its 
Hamiltonian. Thus, in each run the system undergoes 
a unitary evolution, described by a Hamiltonian H^ a '(t) 
which occurs with probability p a . Assuming that the sys- 
tem initially is always prepared in an eigenstate \E a (0)), 
the density matrix for the fluctuating system is given by 
Pit) = Y, a Pa U a{t)Pi Ea ^))Ul(t) . and one finds 



^0 



Tr Tv a B 



/ i 



\E a {0))P\E a (t)) 



Tr 5> Q P ( 



U a \E a (Q))P\E a (t)) 



(18) 

(19) 



For some index a, the application of the AT may fail, 
but averaging over a could mitigate the deleterious ef- 
fects. The exploration of the AT for fluctuating systems 
and mixed states 0] is an important future direction for 
acertaining the validity and limits of the AT. 

In conclusion, we have demonstrated an inconsistency 
implied by the standard statement of the AT and pre- 
sented a counterexample of a two-level system. Both ex- 
amples alert us to the fact that the AT must be applied 
with care. Further work will concern testing the AT for 
various systems, especially those that involve stochastic 
fluctations and mixed states. 

Since this work first appeared as a preprint Q, two 
subsequent preprints appeared that deal with our incon- 
sistency. Sarandy et al. 15] have presented a simplified 



form of the inconsistency which they regard as a valida- 
tion of the standard statement of the AT. We interpret 
their work as an alternative explanation of the cause of 
the inconsistency and a second demonstration that the 
standard statement of the AT, taken as it is, can lead to 
contradictory results. Pati and Rajagopal [lg have found 
a different form of inconsistency associated with the adi- 
abatic GP. Comments on their work and the present in- 
consistency have been made in Ref. [l7j . 
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FIG. 1: Fidelity T between the exact evolution 110H and the 
instantaneous eigenvector for ujq = Is -1 and r = 2-k ■ 10s _1 . 



